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ABSTRACT 

The separation A poak between two peaks in the gamma-ray pulse profile is calculated 
as a function of energy for several polar cap models with curvature-radiation-induced 
cascades. The Monte Carlo results are interpreted with the help of analytical approxi- 
matio ns and discussed in view of the recent data analysis for the Vela pulsar (Kanbach 



1999|) . We find that the behaviour of AP eak as a function of photon energy e depends 
primarily on local values of the magnetic field, i?i OC ai, in the region where electromag- 
netic cascades develop. For low values of Aocai (< 10 12 G), A pcak (e) is kept constant. 
However, for stronger magnetic fields (<; 10 12 G) in the hollow-column model A pcak 
decreases with increasing photon energy at a rate dependent on maximum energy 
of beam particles as well as on viewing geometry. There exists a critical photon en- 
ergy Etum above which the relation A peak (e) changes drastically: for e > Sturm m 
hollow-column models the separation A poak increases (whereas in filled-column model 
it decreases) rapidly with increasing e, at a rate of ~ 0.28 of the total phase per decade 
of photon energy. The existence of critical energy etum is a direct consequence of one- 
photon magnetic absorption effects. In general, etum is located close to the high-energy 
cutoff of the spectrum, thus photon statistics at etum should be very low. That will 
make difficult to verify an existence of etum in real gamma-ray pulsars. Spectral prop- 
erties of the Vela pulsar would favour those models which use low values of magnetic 
field in the emission region (-Bi OC ai I0 11 G) which in turn implies a constant value 
of the predicted A pcak within EGRET range. 

Key words: pulsars: general - pulsars: individual: PSR B0833— 45 - gamma-rays: 
observations - gamma-rays: theory. 



1 INTRODUCTION 

Two prominent peaks are a characteristic feature of gamma- 
ray pulse shapes in the three brightest out of seven gamma- 
ray pulsars detected so far: Crab (PSR B0531+21), Vela 
(PSR B0833-45), and Geminga (J0633+1746). Phase sep- 
aration between the two peaks is very large in each case, 
in the range between 0.4 and 0.5 (e.g. Fierro, Michelson & 
Nolan 1998). The separation, which we denote by A pcak , 
was determined with photons from the entire energy range 
of EGRET. 



Kanbach (1999) suggested that the separation A pcak 



Vela might be energy dependent. The effect would be of the 
order of a few percent or less. The plot of the phase sep- 
aration against energy (fig. 2 of Kanbach 1999) shows that 
A peak decreases by about 5% over 20 energy intervals cover- 
ing the range between ~ 50 MeV and ~ 9 GeV. The scatter 
of points is, however, large enough for this result still to be 
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consistent with the separation staying at a constant level of 
0.43, especially when one rejects two energy intervals: of the 
lowest and the highest value. 

Such effects as suggested by Kanbach can be justified 
qualitatively, at least within polar cap scenarios. Their origin 
may be different at different energy ranges, and their mag- 
nitude may vary as well. For example, Miyazaki & Takahara 
( 1997 ) found dramatic changes in peak-to-peak phase sepa- 
ration due to magnetic absorption effects in their attempts 
to model the Crab pulse shapes. Their numerical calcula- 
tions were performed with low photon energy resolution for 
a model with homogeneous polar cap, and instant accelera- 
tion. 

This new aspect of studying the HE properties of pul- 
sars is potentially attractive. The problem of poor pho- 
ton statistics should become less essential with future high- 
sensitivity missions like GLAST. Then any well established 
empirical relation between the phase separation A peak and 
the photon energy e (including A pcak = const) may help to 
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discriminate in favour of some particular models of pulsar 
activity. 

In this context we present a model of the peak-to-peak 
phase separation in the gamma-ray light curve of Vela and 
confront it with the results of Kanbach (1999). Our aim is 
to present properties of A peak (e) predicted by the polar cap 
model with curvature (CR) and synchrotron (SR) radiation 
being dominant emission mechanisms. This is an exten sion 
of the preliminary results of Dyks, Rudak & Bulik ( 2000). In 



section 2 we outline the model and introduce the input pa- 
rameters for which Monte Carlo simulations were performed. 
Section 3 describes the numerical results and offers their in- 
terpretation; conclusions follow in section 4. 



2 THE MODEL 

The presence of two peaks with large (0.4 - 0.5) phase sep- 
aration in gamma-ray lightcurves within single polar cap 
models requires a nearly aligned rotator (e.g. Daugherty & 
Harding 1994) where the following three characteristic an- 
gles are to be of the same order: a - the angle between the 
spin axis ft and the magnetic moment fj,, # 7 - the open- 
ing angle between the direction of the gamma-ray emission 
and fj,, and £ - the angle between ft and the line of sight. 
For a canonical polar cap and instant electron acceleration, 
# 7 roughly equals 0.02/\/P radians only (where P denotes 
a spin period). To avoid uncom fortab ly small characteristic 
angles, Daugherty & Harding ( 1996) ) , postulated that pri- 
mary electrons come from extended polar caps, and with 
the acceleration occuring at a height h of several neutron- 
star radii R us . The latter assumptio n ma y be supported by 
the results of Harding & Muslimov ( 1998] ) who investigated 
in a self-consistent way particle accel eration by th e electro- 
static field due to field-line curvature (|Arons 1983 Jand iner- 
tial frame dragging effect ( Muslimov fc Tsygan 1992j ). Hard- 
ing & Muslimov found that a stable accelerator, with dou- 
ble pair-formation-front controlled by curvature radiation, 
is possible at a height of about 0.5 to 1 stellar radii. Here 
we use a polar cap model combined with the assumption of 
a nearly aligned rotator. Most ingr edien ts of the pc-model 
come from Daugherty & Harding ( 1982 ). Geometry of the 
magnetic field of a neutron star is assumed to be well de- 
scribed by a static, axisymmetric dipole. 

Within a given polar cap (pc) model with a fixed value 
of a there are two possible values of viewing angle £ re- 
sulting in an identical peak separation A peak (defined as a 
fraction of 2n) but with a reversed order of the two peaks 
(in terms of a leading peak, and then a trailing peak). Fig. [j] 
illustrates the ambiguity in the definition of a double-peak 
pulse. A simple geometrical relation connects the three an- 
gles of interest - a, £ and # 7 - to the phase separation 

^peak . 



COS 



cos S sin a sin £ + cos a cos £, 



(1) 



(eg. Ochelkov & Usov 1980), where S = 7rA pcak for Ciargc, 
and 5 = 7r(l — A pcak ) for (small This relation holds as far 
as aberration of photon propagation due to rapid rotation 
is neglected. 

Throughout the paper we always take the case of the 
larger £ in each model (4 models are considered). Its value 
(Ciarge = 3.75,4.5,15., and 3.65 degrees for models A, B, 
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Figure 1. The double-peak pulses at 100 MeV (thin line) and 500 
MeV (thick line) calculated for model A. Depending on viewing 
angle £ between the spin axis f2 and the line of sight each single 
polar cap model yields two types of double-peak pulses. The upper 
and lower panels show the pulses for two values of the angle (: 
Csmall = 2.5 deg and f largc = 3.75 deg, respectively. 



C, and D, respectively) along with our choice for the angle 
a (see Table 1) is to yield the separation A peak = 0.43 at 
300 MeV. 

The dominant HE emission process is the curvature ra- 
diation (CR) by ultrarelativistic beam particles (primary 
electrons which leave the polar cap), followed by magnetic 
pair production with the subsequent synchrotron emission 
(SR). Our numerical code to follow the cascades induced 
by beam particles is based on Daugherty 



Harding (|982J), 
and takes advantage of the following approximations rele- 
vant to the problem of directional and spectral distributions 
of the photons: The curvature photons are emitted tangen- 
tially to the local magnetic field direction in a frame rotating 
with the star. The created e ± -pairs are assumed to follow 
the directions of their parent CR photons and they share 
the photons' energy equally (for justification see Daugh- 
erty & Harding 1983). The synchrotron photons are emitted 
perpendicularly to the local magnetic field direction in a 
frame comoving with electron/positron center of gyration. 
The pairs are assumed not to change their position on the 
field line when radiating (for magnetic field strengths con- 
sidered here, energy-loss length scales due to SR are of the 
order of 10~ 6 cm). The emergent high-energy spectrum is 
a sup erposition of CR and SR. We follow Rudak & Dyks 
to calculate detailed broad band energy spectra of 



( [1999D 



the high-energy emission. 

Beam particles are injected along magnetic field lines 
into the magnetosphere either from the outer rim of the 
polar cap (hollow cone column) or from the entire polar 
cap surface (filled column). We use two simple models for 
their acceleration to ultrarelativistic energies: 1) instant ac- 
celeration and 2) acceleration due to a uniform longitu- 
dinal electric field £ over a scale height Ah. The pulse 
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Table 1. Model parameters. 







a 




#init 












primary electrons 




[10 12 G] 


[deg] 


[Rns] 


[0pc] 


Model A 


1.0 


3.0 


0.0 


1.0 


E; n i t = 8.68 TeV, no acceleration 


Model B 


1.0 


5.0 


1.0 


1.0 


E ln it = 20.0 TeV, no acceleration 


Model C 


3.0 


10.0 


2.0 


2.0 


^init = 0.5 MeV, acceleration (see eq.(^|)) 












with V = 1.5 X 10 13 Volts 


Model D 


1.0 


3.0 


0.0 


[0, 1] 


E; n i t = 8.68 TeV, no acceleration 



shapes as a function of photon energy were calculated nu- 
merically for four sets of initial parameters (hereafter called 
models A, B, C and D). Table 1 features most important 
model parameters. In models A and B the primary elec- 
trons are distributed evenly along a hollow cone formed by 
the magnetic field lines from the outer rim of a canonical 
polar cap, i.e. with a magnetic colatitude ln i t ~ pc , where 
6 pc ~ (2-7T R ns /c P) 1 / 2 radians at the stellar surface level 
(h = 0). The beam particles are injected at a height /linit 
with some initial ultrarelativistic energy E- m it (listed in Ta- 
ble 1) and no subsequent acceleration takes place. The main 
difference between models A and B is due to different values 
of /linit (equal to and 1 R ns respectively) which result in 
different locations of origin of secondary particles. Changing 
these locations is an easy way to modify spectral properties 
of emergent radiation and enables to change (preferrably - to 
increase) the angle a as constrained by t he observed ene rgy- 
averag ed peak separation A peak w 0.43 ( [Kanbach 1999 ). In 
model C we assume ft n it = 29 pc . The primaries are injected 
at /linn = 2i? ns with .Emit = rnc 2 and undergo acceleration 
by a longitudinal electric field £ present over a characteris- 
tic scale height Ah = 0.6 R ns , resulting in a total potential 
drop Vo: 



£(h) 



f V /Ah, for /linit <h< (/linit + Ah) 

\ 0, elsewhere. 



(2) 



For comparison, we considered model D with a uniform 
electron distribution over the entire polar cap surface (i.e. 
#init £ [0, #p C ]). All its remaining features are indentical to 
model A. 

The values of .Emit in models A, B and D, and the poten- 
tial drop Vo in model C were chosen to yield similar number 
of secondary pairs — about 10 3 per beam particle. In all 
cases the spin period of the Vela pulsar P = 0.089 s was 
assumed. 



3 RESULTS 

General properties of the function A peak (e) may be de- 
scribed in short in the following way: In the low-energy range 
(i.e. below a few GeV) the peak separation either remains 
constant with the increasing photon energy (Models B, C, 
and D) or slightly decreases (Model A). In either case, how- 
ever, the slope of A pcak versus e looks consistent with the 



results of Kanbach (1999). Then, around a few GeV there 
exists a critical energy etum, at which the separation A peak 
undergoes a sudden turn: in the high-energy domain (i.e. for 
e > £tm-n) it increases in our hollow-column models (A, B, 
and C), whereas it decreases in the filled-column model (D), 
at a rate ~ 0.28 phase per decade of photon energy. The ex- 



istence of £tum is a direct consequence of magnetic absorp- 
tion (7B — > e*) in the magnetosphere (see also Miyazaki 
& Takahara 1997). [Note, that etum is not equivalent to a 
high-energy cut-off in a spectrum due to the magnetic ab- 
sorption] . 

The dependence of peak separation upon photon energy for 
all four models is presented in Fig.^ (upper panels). 

To understand why the slope in the low-energy domain 
(e < £tum) is different in different models, we will now 
present several factors which are of interest in this respect. 

Let us start with possible consequences of the orien- 
tation on A pcak (e) in this context, i.e. with the choice of 
the angles a and £. For a given emission pattern in a 
frame of magnetic dipole rotating with the star and for 
a set of values of angles a and £ fulfilling the condition 
A peak (300 MeV) = 0.43, a line of view can cross the hol- 
low cone of emission at different impact angles. Let us esti- 
mate how geometry alone would affect the relations A pe (e) 
shown in Fig. |^. We neglect aberration of photon propa- 
gation due to rotation when calculating directional emis- 
sion pattern as seen from some inertial observer frame. For 
A pcak staying close to 0.5 (roughly between ~ 0.3 and ~ 0.7) 
cos(7rA pcak ) ~ J - 7rA pcak , the value of costf 7 is approxi- 
mately proportional to A poak with a constant of proportion- 
ality equal to — 7rsinasin£. Moreover, for a radiating parti- 
cle sliding along a dipolar field line with the dipole constant 
k = r smT 2 8 (r, 8 are the coordinates of the particles in 
the dipolar frame), one can link the opening angle # 7 with 
a radial position r of the particle: 



cos 



V^ 



(3) 



Since in polar cap models ^ < 1, the right-hand side of 
eq. ([]) can be approximated by cos # 7 ~ 1 — | j£ which along 
with eq.dlt) gives 



8 

r ~ — . 
9 

where 



TvaA 



peak 



+ 6 



(4) 



k sin a sin £ , 



and 
6 



^ k ( 1 — — sin a sin ( — cos a cos () . 



9 v 2 

The slope of the relation A pcak (e) is a combination of two 
factors: 1) viewing geometry and 2) directional and spectral 
changes in the radiation yielding the peaks with increasing 
radial coordinate r: 



dA pcak _ dA pcak dr 
de dr de 



(5) 
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Figure 2. Upper panels: Phase separation A pcak versus photon energy e of the emission peaks found with Monte Carlo calculations 
for models A, B, C and D (dots). The overall observational trend, taking into account a substantial scatter of points (see fig. 2, middle 
panel, of Kanbach 1999), is marked schematically for reference in all four panels as a parallelogram. The assumed values of £ - the angle 
between CI and the line of sight - are also indicated. The steep dashed lines in the three left-hand side panels (hollow-column models) 
are the approximate analytical solutions for A poak determined by magnetic absorption effects (7B — > e^). They reproduce the Monte 
Carlo results with good accuracy in the high energy limit (the GeV domain). 

Lower panels: Energy output per logarithmic energy bandwidth at the first peak as a function of photon energy e for models A, B, C 
and D. The synchrotron (long dashed) and the curvature (short dashed) components are marked for models A, B, and C. [Note: these 
are not instantaneous spectra.] Units on the vertical axis are arbitrary. 



Any changes in viewing geometry (i.e. the angles a and £) 
will affect the slope dA pcak /cfe via <fA pcak /dr = the 
latter being calculated from eq.(^). For example, for models 
A and B the parameter a equals 1.46i? ns and 2.91J? ns , re- 
spectively. Should dr/de be identical for these two models, 
the slope in model A would be about two times bigger than 
in model B. Fig. | (upper panels) does not show any such 
effects: the slope in model B is ~ for e < etum- This is 
because the viewing geometry effects are dwarfed by differ- 
ences in directional and spectral properties of the radiation 
in these models. 

To analyse directional and spectral properties of the 
peak emission it is instructive to begin with a simplified, one- 
component model. Suppose that the only contribution to the 
emission is due to optically thin CR, i.e. both magnetic ab- 
sorption and synchrotron emission by secondary pairs cre- 
ated due to this absorption are neglected. Numerical simula- 
tions of this case (not included in this paper) show that the 
peak separation in such a case, which we denote as A pcak , 
stays constant as a function of e for the models with instant 
acceleration (A, B, D), while it increases with increasing e 
for model C (acceleration over Ah) at a rate dependent on 
the electron acceleration rate. The former case can be un- 
derstood by arguing that energy E = jmc 2 of the electron 
decreases monotonically from its starting value -Emit, while 
its radius of curvature p cr increases. Recalling the properties 



of the curvature continuous spectrum due to a single particle 
(e.g. Landau & Lifshitz 1973) this makes the contribution 
to the spectrum per unit distance by the electron to be the 
highest one just at the initial altitude /i m it- In consequence, 
the contribution to the spectrum per unit phase angle by a 
bundle of electrons moving along a set of open field lines is 
the highest one at ftinit. Therefore, it is the opening angle 
of the CR photons at /linit which fixes the phase separa- 
tion A P ° ak of the two peaks regardless the photon energy. 
In the latter case the qualitative behaviour of A P ° ak (e) can 
be understood by assuming that the phase of the pulse peak 
at a given photon energy e is approximately determined by 
such altitude at which accelerated electrons reach the en- 
ergy ■ymc 2 which satisfies the condition 0.29 £ cr (7, Per) ~ e, 
where e cr = |c?i7 3 p^. 1 is a characteristic energy of the CR 
spectrum and p cr is a local radius of curvature; both 7 and 
p cr are functions of altitude h. 

We now relax these two simplifications and present the 
consequences. First, the inclusion of magnetic absorption re- 
sults in a strong response of A P r cak at highest photon ener- 
gies - between some critical energy etum and the high-energy 
cutoff, as mentioned in the first paragraph. This effect will 
be addressed in the final part of this section. Second, adding 
up synchrotron component (SR) due to pairs changes the 
properties of high-energy emission for s < etum (i-e. below a 
few GeV) significantly. Most notable is the domination of SR 
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Figure 3. Phase separation A peak versus photon energy e corresponding to the directions of propagation of curvature photons emitted 
by electrons located on annuli (each annulus corresponds to a given radial coordinate and a dipolar constant; it translates into A peak for 
a fixed set of angles between the spin axis n, the magnetic moment fi, and the line of sight). The curvature photons are here assumed 
to have a characteristic energy e cr = 1 .bch-y^p^. 1 (the dimensionless energy 7 e ] of the electrons as well as the local radius of curvature 
p cr change along electron's trajectory in a model dependent way). Dashed lines (as in Fig. ^. upper panels) represent the limiting energy 
£ csc of curvature photons for which the magnetospere becomes opaque. Photons located to the right of the dashed line in each panel 
convert into eP 1 pairs. This part of A peak vs. e is shown as dotted line. The pairs produce the synchrotron spectral component which in 
turn determines the character of A pcak vs. e in the low-energy domain (see Fig. El upper panels). Note a wide range of A pcak values for 
which curvature photons are absorbed in Model A, in contrast to Models B and C. 



over CR in terms of intensity. The total energy output per 
logarithmic energy bandwidth at the first peak as a function 
of photon energy e for models A, B, C and D is presented 
in Fig. |^ (lower panels). Both components - SR and CR - 
are marked schematically to show their relative importance. 
Consequently, A pcak is affected by directional and spectral 
properties of the SR emission. 

The behaviour of A pcak is shown in Fig. [| (upper pan- 
els). Its slight decrease with increasing e (model A) or no 
change at all (models B, C, and D) - is due to a combina- 
tion of factors which determine the directional and spectral 
properties of SR. These include energy and pitch angle dis- 
tributions of secondary e pairs, as well as their vertical 
spread within the magnetosphere combined with a strength 
of the local magnetic field. These factors change from one 
model to another: 

Model A - In strong local magnetic fields (Biocai ~ 10 12 
G), efficient pair production requires lower electron ener- 
gies than in a low- field case (like in model B). The produc- 
tion occurs over a wider range of altitudes. Consequently, 
the spectrum of SR contributed locally by the pairs evolves 
considerably over this range of altitudes. Lorentz factors of 
gyration j± are of low values (of the order of mc 2 /(15eb), 
where eb = mc 2 B\ oca i / B cr it) • The resulting local SR spec- 
tra are, therefore, very narrow (see also the lower panel of 
fig. 1 in Rudak & Dyks 1999): they spread between the (lo- 
cal) characteristic SR energy e sr and the (local) cyclotron 
turnover energy e ct = 1.5es/ sin ip (where ip is the pitch an- 
gle) which roughly satisfy e sr /e c t ~ 7x, and the ratio does 
not exceed 10 3 in model A. With increasing height (there- 
fore, with decreasing Bi oca i) both e sr and e c t move towards 
lower and lower values. The final effect of this softening is 
then a notable decrease of A pcak with increasing photon 
energy (Fig. |[ upper panel, Model A). In even stronger 
local magnetic fields (but not exceeding B cr it) the rate of 
decrease of A peak with increasing e may be much larger, 



because e -pairs are produced with extremely low j± and 
synchrotron/cyclotron photons emitted at any altitude con- 
centrate near the local cyclotron energy eb- 

Model B - When cascades are to develop in a relatively weak 
magnetic field, -Bi oca i 10 11 G, very high electron energy 
Sinit is required. The high value of .Emit means rapid CR 
cooling which brings the electron energy quickly below the 
level required for pair production. Curvature photons be- 
come too soft for pair creation via magnetic absorption very 
quickly after injection. The bulk of pairs is created by elec- 
trons with their energies confined to a narrow range at E- m it, 
and consequently the SR component originates within a nar- 
row range of magnetospheric altitudes (radial positions r). 
The resulting A pcak does not change with photon energy e. 

The difference between model A and B in radial exten- 
sion r of the regions of origin of curvature photons which are 
absorbed, producing pairs and eventually SR, is presented in 
Pig. |^ (right-hand vertical axes). The difference is more ap- 
pealing when presented in terms of A peak (left-hand vertical 
axes) . 

Model C - After initial stage of linear acceleration, the elec- 
trons enter a regime of 'radiation reaction limited accelera- 
tion'. Over a considerable range of altitudes (3.3i? ns Ss r Si 
3.6-Rns) the electrons' energy remains approximately con- 
stant and so does the pair production efficiency (it is equal 
to ~ 4 • 10 pairs per centimetre of the primary electron 
path). In such conditions, any evolution of the SR-spectrum 
over this range might affect A peak as a function of e. Never- 
theless, no significant changes occur in A pcak (e). This is be- 
cause a single-particle SR-spectrum almost does not change 
its shape within the range of altitudes with stable, strong 
pair production. The balance between the acceleration rate 
and the CR-cooling rate stabilizes e cr as well as e sr . Any 
spectral changes of SR in its low-energy part, nearby e c t, are 
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not relevant in the context of gamma-rays, since they occur 
near 100 keV, i.e. well below the energy range of EGRET. [|] 
Model D (a uniform distribution of primary electrons over 
the polar cap but otherwise identical to model A) - the final 
effect does not resemble the decrease of A pea with increas- 
ing photon energy found for model A. Apart from numerical 
fluctuations (see Fig. ^, upper panel) the separation A pcak 
remains approximately constant as a function of e. 

Finally, let us discuss the inclusion of magnetospheric 
opacity due to 7B — > e*. This effect becomes important 
above ~ 1 GeV, where most power is due to CR in all mod- 
els. In this regime the position of each peak in the pulse 
is determined by magnetic absorption, and this results in a 
strong response of A pcak between eturn ii 1 GeV and a high- 
energy cutoff: At e = £tum the separation A peak undergoes 
a sudden turn and starts increasing (or decreasing) rapidly 
for e > eturn- For the hollow-column models (A, B and C) 
the photons in both peaks of a pulse come from low magne- 
tospheric altitudes with narrow opening angles. When e is 
high enough these photons will be absorbed. Photons which 
now found themselves in the 'new' peaks come from higher 
altitudes (the magnetosphere is transparent to them) and 
have wider opening angles. In other words, inner parts of the 
'original' peaks in the pulse will be eaten-up and the gap be- 
tween the peaks, i.e. the peak separation A pcak , will increase 
(FigH). Our Monte Carlo results for A pcak at e > £, nm may 
be reproduced with good accuracy by a simple analytical 
solution of t 1 b > 1 which has to be combined with eq.(Q). 
The requirement r 7 _g > 1 is particularly simple (with some 
well known approximations being used) since it refers to a 
photon created with a momentum tangential to local dipo- 
lar magnetic field line at a height h above the neutron star 
surface (at radial coordinate r = h + R ns )- The photon will 
undergo magnetic absorption (to be more precise, with a 
probability of [1 — exp(— t 7 b)]) if its energy e satisfies the 
following condition: 



e > £ c 



(6) 



(cf. eq. 11 in Bulik et al. 1999). This formula is valid for 
hollow-column models but may be used for any dipolar field 
line (the factor P 1 / 2 comes just from choosing the outer rim 
of a polar cap to be the site of field- line footpoints). Al- 
though the pulsar spin (which leads to an aberration and 
a slippage of magnetic field under the photon's path) has 
been neglected in derivation of eq. ^ , the formula gives e osc 
in excellent agreement with Monte Carlo method for the 
emission regions placed up to several R na above the surface 
and rotation periods typical for strong-field pulsars (<; 10~ 2 



T The behaviour of SR in its low-energy part (i.e. in hard X-rays) 
may then lead to increase of A peak in the hard X-ray domain in 
comparison to A pe in the gamma-rays. Assuming alternative 
to our choice definition of A pcak , with Csmall (see Section 2), the 
behaviour of A pcak would be a mirror-reflection of our case. Then 
this effect might explain low value of A peak withi n the 2 — 30 keV 
range found by Strickman, Harding & dejager (1975) in RXTE 



data of Vela. Moreover, below 2 keV the synchrotron component 
in the model spectrum drops below the level of CR-emission, and 
the value of A peak as for the gamma-rays is expected to be re- 
sumed there. 



s). One may look at the eq.(^|) as the condition for a ra- 
dius of escape r csc at a given energy e. Photons of energy 
e will escape the magnetosphere only if they are emitted at 
r > r csc which satisfies e OS c(r) > e. The radial coordinate 
r e sc(e) has to be combined now with eq.(^) to give A poak rel- 
evant for the 'magnetic absorption 'regime. This analytical 
A pcak is shown as dashed lines in Fig^ (upper panels) and 
in Fig.^| whereas the filled dots are the Monte Carlo results. 
This branch of solution intersects the horizontal line set by 
A poak = 0.43 at etum ^ 0.9, 4.5, and 3 GeV for models A, 
B, and C, respectively. 

For model D, with a uniform distribution of primary 
electrons over the polar cap (but otherwise identical to 
model A) the changes of A pcak above eturn, occur in the 
opposite sense. Unlike in previous models, here both peaks 
of the pulse are formed by photons emitted tangentially to 
magnetic field lines attached to the polar cap at some open- 
ing angle 9- m it < 6 pc ■ These photons are less attenuated than 
those coming from a hollow-column, and in consequence the 
peak separation drops. A sim ilar behaviour was obtained by 
Miyazaki & Takahara (1997), in their model of a homoge- 
neous polar-cap. 

We have also investigated the behaviour od A pea (e) 
above eturn for other distributions of primary electrons over 
the polar cap. We have considered intermediate cases be- 
tween models A and D, (i.e. with a uniform coverage of only 
an outer part of the polar cap area between some inner 
radius r in < r pc and the polar cap radius r pc ), and mod- 
els with uniformly filled interior of the polar cap surface 
but with increased electron density along the polar cap rim 
(cf. Daugherty & Harding 1996). We conclude that regard- 
less the actual shape of the active part (i.e. 'covered' with 
primary electrons) of the polar cap (either an outer rim, or 
an entire cap, or a ring, or entire cap/ring + increased rim 
density), one does expect in general strong changes in the 
peak separation to occur at photon energies close to high- 
energy spectral cutoff due to magnetic absorption. 

A word of technical comment seems to be appropriate 
for A pcak = 0.43 at ~ 3 00 M eV. It appears that a technique 
adopted by Kanbach (1999) of fitting the observed pulse 
shapes with asymetric Lorentz profiles tends to overestimate 
the true value of A pea by a few thousandth parts of phase. 
Therefore the actual value may be closer to 0.42 than 0.43 
(Maurice Gros, private communication). Nonetheless, this 
shift in A peak does not change any conclusions of this work. 



4 SUMMARY 

In thi s pap er we addressed a recent suggestion of Kan- 
bach (11993) that peak separation A peak in the double-peak 
gamma-ray pulses of the Vela pulsar may monotonically de- 
crease with increasing photon energy at a rate ~ 0.025 phase 
per decade in energy over the range 50 MeV to 9 GeV, We 
calculated gamma-ray pulses expected in polar-cap models 
with magnetospheric activity induced by curvature radiation 
of beam particles. Two types of geometry of magnetospheric 
column above the polar cap were assumed: a hollow-column 
associated with an outer rim of the polar cap and a filled 
column associated with a uniform polar cap. Four models 
were considered with two scenarios for the acceleration of 
beam particles. Pulsed emission in the models was a super- 
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position of curvature radiation due to beam particles and 
synchrotron radiation due to secondary pairs in magne- 
tospheric cascades. The changes in the peak separation were 
investigated with Monte Carlo numerical simulations. 

We found that regardless the differences in the mod- 
els, the peak separation A pcak below a few GeV, where the 
emission is dominated by synchrotron component, is either 
a weak decreasing function of photon energy e, or remains 
constant. Both variants may be co nsider ed to be in agree- 
ment with the results of Kanbach (1999) for the latter are 
affected by large statistical errors. A particular behaviour of 
A peak depends on a combination of several factors, includ- 
ing strength of magnetic field in the region of pair formation 
and model of electron acceleration (both of which determine 
spectral and directional properties of the radiation at differ- 
ent altitudes), as well as viewing geometry. Essentially, in 
strong fields, Bi oca i <; 10 12 G, A peak decreases with increas- 
ing photon energy e, whereas for Bi oca i < 10 12 G, the peak 
separation A pcak stays at a constant level. 

Moreover, we found that due to the magnetic absorp- 
tion (7B > e*) there exists a critical energy E lurn at which 
the peak separation A peak makes an abrupt turn and then 
changes dramatically for e > etum. It increases in the hollow- 
column models (A, B, and C) and decreases in the filled- 
column model (D), at a rate ~ 0.28 phase per decade of pho- 
ton energy. The numerical behaviour of A peak in this regime 
in the hollow-column models was easily reproduced to high 
accuracy with a simple analytical model of magnetospheric 
transparency for a photon of energy e, and its momentum 
tangential to local dipolar magnetic field line at a site of its 
origin. An exact value of etum is model-dependent but it is 
confined to a range between ~ 0.9 GeV and ~ 4.5 GeV. 

To find such a hypothetical turnover of A pcak in real 
observational data would require, however, high-sensitivity 
detectors, since for e > etum the expected flux of gamma- 
rays drops significantly. If detected, this turnover would be 
an important signature of polar cap activity in gamma-ray 
pulsars. It would support the notion that high-energy cut- 
offs in gamma-ray spectra of pulsars are due to magnetic 
absorption. 

The CR-induced cascades models, like those considered 
in this work, are not the only possibility for nearly aligned 
rotators to produce double-peak pulses with large phase sep- 
arations. There exists an alternative class of models - with 
pair cascades above polar cap induced by magnetic inverse 
Compton scatterings (ICS) of primary electrons in the field 
of soft photons from the stellar surface - proposed in a se- 
ries of papers (e.g. Sturner & Dermer 1994, Sturner et al. 
1995). In particular, Sturner et al. (1995) present the de- 
tailed Monte Carlo model spectra of the Vela pulsar. They 
also present pulse profiles at a fixed energy of 50 MeV (for 
several viewing angles) but no word of comment is given re- 
garding the problem of A pcak versus photon energy. We ex- 
pect the outcome to be qualitatively similar to our results. 
First, the scatterings take place mostly within a very limited 
height above the polar cap surface (below h ~ R P c) and the 
preferred directions of propagation of the ICS photons will 
be fixed by magnetic field lines just above the surface. There- 
fore, A peak due solely to ICS photons should stay constant 
for a wide range of energy. Inclusion of synchrotron pho- 
tons due to pairs is unlikely to notably affect A poak unless 
the pair formation front is vertically more extended than for 



CR-induced cascades. Second, some turnover point at etum 
not exceeding 1 GeV should be present due to magnetic 
absorption. The behaviour of A pcak for e > etum should 
roughly follow the dashed lines in Fig. 2 (upper panel) and 
Fig. 3 as long as the assumption about photons (which are 
to be absorbed) propagating tangentially to local dipolar 
magnetic field line at their site of origin remains valid for 
majority of ICS photons. 

To verify this qualitative picture would, however, require 
detailed numerical calculations. 
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